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We study how the dynamical vector degree of freedom in modified gravity affects the CMB B- 
mode polarization in terms of the Einstein-aether theory. In this theory, vector perturbations can 
be generated from inflation, which can grow on superhorizon scales in the subsequent epochs and 
thereby leaves imprints on the CMB B-mode polarization. We derive the linear perturbation equa- 
tions in a covariant formalism, and compute the CMB B-mode polarization using the CAMB code 
modified so as to incorporate the effect of the aether vector field. We find that the amplitude of the 
B-mode signal from the aether field can surpass the contribution from the inflationary gravitational 
waves for a viable range of model parameters. We also give an analytic argument explaining the 
shape of the spectrum based on the tight coupling approximation. 
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I. INTRODUCTION 

It is widely believed that the present Universe is dom- 
inated by the two dark components: cold dark matter 
and dark energy. Dark matter plays an essential role in 
explaining galaxy rotation curves and in structure forma- 
tion, while dark energy is presumably responsible for the 
current cosmic acceleration. The presence of the dark 
components is thus perceived through the gravitational 
interaction, having not yet identified what they really 
are. It is therefore legitimate to think of these major 
mysteries of today's cosmology as a mystery of gravity. 
This motivates us to explore long-distance modification 
of the gravitational law, asking to what extent general 
relativity (GR) is correct on cosmological scales. 

Modification of gravity is most commonly made by 
adding an extra scalar degree of freedom a la Brans- 
Dicke gravity In recent years, various refined models 
of scalar-tensor gravity have been proposed which would 
be an alternative to material dark energy while being 
consistent with solar-system constraints. They include 
chameleon /(i?) [l, Q and Galileon 0, Q theories, and 
have been tested against cosmological observations [y| . It 
is also possible to modify the spin-2 sector as in massive 
gravity and bi-gravity theories Q. 

In this paper, we are going to consider a hypothetical 
vector degree of freedom of gravity. Specifically, we shall 
focus on the Einstein-aether (EA) theory proposed by Ja- 
cobson and Mattingly Q, in which a fixed norm vector 
field with a Lorentz-violating vacuum expectation value 
takes part in the gravitational interaction. The effect of 
the aether on the cosmological background was clarified 
in [l^ . The scalar cosmological perturbations in the EA 
theory and their impact on the cosmic microwave back- 
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ground (CMB) temperature anisotropy have been stud- 
ied in [nMl^ . Recently, Armendariz-Picon et al. [l5| 
performed a comprehensive analysis on cosmological per- 
turbations in the EA theory. See also Refs. for 
the other aspects of the EA theory, such as spherically 
symmetric solutions and compact objects. Interestingly, 
it was recently pointed out that the healthy extension [2^ 
of Horava's quantum theory of gravity [2l| reduces to a 
special case of the EA theory at low energies [l^l • Cosmo- 
logical perturbations in the healthy extension of Horava 
gravity were studied in [23| ■ 

The purpose of the present paper is to clarify the im- 
pact of the aether vector field on the CMB polariza- 
tion. The CMB polarization arises from all the three 
types of cosmological perturbations, i.e., scalar, vector, 
and tensor perturbations. Among them, as pointed out 
by the vector perturbations most effectively gener- 
ate the B-mode polarization. However, the effect of vec- 
tor perturbations has been less investigated because the 
vector mode decays unless sourced, 6.17., by topological 
defects [2^, [2^ , while the scalar and tensor modes are 
certainly generated from inflation Other possible 

ways of seeding vector perturbations include the neutrino 
anisotropic stress [2^, the second-order effect [2£' 
and primordial magnetic fields generated somehow 
Modifying the vector sector of gravity offers a yet an- 
other possibility of producing vector perturbations, leav- 
ing a unique signature in the CMB polarization due to 
nontrivial dynamics of the aether field. 

The paper is organized as follows. In the next section 
we introduce the EA theory and the basic equations. In 
Sec. Ill, we describe the dynamics of vector perturbation 
in the EA theory using the covariant approach, in or- 
der to incorporate the aether vector field into the CAMB 
code. We then specify the initial conditions for the per- 
turbation evolution in Sec. IV. Our numerical results are 
presented in Sec. V. In Sec. VI, we examine the spectrum 
shape in an analytic approach using the tight coupling 
approximation. We draw our conclusions in Sec. VII. We 
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will use the sign convention (- 



III. COVARIANT APPROACH 



II. EINSTEIN-AETHER THEORY 

The action of the EA theory is given by 
A'P r 

5 = ^ y d^xV^in + Ca] + S,n, (1) 

where 

Ca = -[ciV,AVMb + c2(VbA^)2 + c3Va^''V6^" 
+C4A''A''VaA''VbA^] + X{AbA'' - 1) (2) 

is the Lagrangian of the aether field and Sm is the action 
of ordinary matter. It is assumed that the aether is not 
coupled to the matter field directly. 

Variation with respect to the metric leads to the Ein- 
stein equations 
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where TZab is the Ricci tensor, n = Mp^, 



Tab -^^Agab 



6C, 

5g' 



ab 



(3) 



(4) 



is the energy-momentum tensor of the aether, and Tab is 
the energy momentum tensor of ordinary matter. Ex- 
plicitly, we have 

Tab = Ci[{yaA'){VbAc) - {VcAa){V''At)] 

-HCiVe [A'V^aA^ + (V^^(J^) - ^(aVb)A^] 
+C2ffabVc {A''\/dA'^) 

-CiA'=A''{WcAa){WdAb) 

-C4V, [AaAtA^VaA^ ~ 2A- A\V ^A^^A^)] 



--^^Agab + >^AaAb. 



(5) 



Ordinary matter includes photons, baryons, etc., so that 
we write Tab = '''ab ' where i labels different compo- 
nents. Note that Ca is taken to be a functional of A°' 
rather than Aa- Variation with respect to A" yields the 
equation of motion for the aether: 

+C4 [Vb {A'^A'VcAa) - A\Vi,A,)\7aA'=] = -AA,. 

(6) 

Finally, variation with respect to the Lagrange multiplier 
A gives the fixed norm constraint 

AaA- = 1. (7) 

It is convenient to use the following abbreviations: 

Cl3=Ci+C3, Ci4=Ci+Ci, 

a = Ci+3c2+C3, Ci23 = Ci + C2 + C3 (8) 



To describe background cosmology and the evolution of 
vector perturbations, we employ the covariant equations 
obtained by the method of 3 -I- 1 decomposition. We 
begin with splitting physical quantities with respect to 
observer's 4- velocity u"". Following the usual procedure, 
the projection tensor is defined as hab '■= gab — UaUb and 
the covariant spatial derivative acting on a tensor 
field Tj!::: is defined as D'^T^'';;; := hfh''j • • • /ij • • • V'T^.;;. 
The energy-momentum tensors of ordinary matter and 
the aether field are decomposed respectively as 



p^'^UaUb-p^'^hab + 2q\yb) + Kt (9) 

(10) 



Tab = PUaUb - phab + '2q{a'Ub) + TTab, 

while VaWf) is decomposed as 



VaUb = -9hab + (Tab + i^ab - Ua.Ub- 



(11) 



Here, aab ■— ^{a^b) ~ (l/3)VcW'^/iah is the shear tensor, 
LOab ■= D[aW6] is the vorticity, 6 := VqU° is the expansion, 
and the overdot denotes time derivative ' := u^Va. The 
expansion may be written as 9 = 3S/S, where S is the 
averaged scale factor. 

The conservation equations for the matter energy- 
momentum tensor imply 



p + d{p + p)+ D"qa = 0, 



(12) 



Qa + ^Oqa + ip+p)Ua- ^aP + ^''iTab = 0, (13) 



where p J2iP^'\ P la = T,t1a \ and nab = 

Ei'^'ab- ^^^'^ have the corresponding conservation 
equations for the aether: 



^ + e{p+p)+D''qa=0, 



(14) 



L + ^Oqa + {p+ p)Ua - ^aP + ^^^ab = 0. (15) 

Let us first consider the homogeneous and isotropic 
limit. In this limit, the constraint ([7]) implies A° = u°- . 
Substituting this to the equation of motion ([6]) , one finds 



(16) 



The energy density and pressure of the aether are then 
given by 



P 



C2 



A - 



P 



--{29 ' ^2 
6 



(17) 
(18) 



It can be seen that p and p are expressed in terms of the 
expansion 6, and they are the same as what appear in 
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the left hand side of the Einstein equations. This means 
that the effect of the aether on the background evolution 
is just to renormalize the gravitational constant: k — )• 
K := (1 — a/2)~^K [l^. The background equations are 
thus given by 



n' = --s\p + 3p), 

6 



(19) 
(20) 



where we have introduced the comoving Hubble param- 
eter, 7i := S6/3, and the derivative with respect to the 
conformal time, ' := Su°'Wa. 

In order for the background Friedmann equation (jl9p 
to have a solution, the condition 



a < 2 



(21) 



must be satisfied for positive matter energy density and 
positive K. On the other hand, the effective gravitational 
constant on small scales kat is also different from the bare 
one: kat = (1 + ci4/2)-1k [13. The difference between 
these two effective gravitational constants is constrained 
by nucleosynthesis as |1 — k/kjv| < 10%. In terms of the 
aether parameters, this constraint is roughly expressed 



ci4 + a< 0.2. 



(22) 



Note that this is trivially satisfied if we consider the spe- 
cial case a = —cn. Actually, this is the special combi- 
nation for evading the existing observational constraints 
on Ci (see Appendix A), and we will often use this case 
later. 

Having studied the background effect of the aether, we 
then move on to the dynamics of vector perturbations. 
As employed in [28j . we choose Ua to be hypersurface 
orthogonal, so that curl Ub = Uft = at linear 

order. This simplifies the following analysis. 

At linear order, the aether field can be written as 



Ah 



(23) 



where y^*-* and Vb are first order quantities. V^'^'^ cor- 
responds to a scalar perturbation which we do not con- 
sider in this paper, while Vb a vector perturbation that 
satisfies BbV^ = 0. The fixed norm constraint ([7]) 
leads to UbV^ = 0, i.e., V'' is a spatial vector. Since 

VaVh DaVb - il/3)eVaUb + UaVb, WC havC 
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tOK, 



UVaUb + UaVb. (24) 



Substituting this to Eq. 
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^aVb 

)]), we find, up to first order 



Cl3 



+ 0' K 



- ClDVa 



\Ua + \Va, 



C29Ua 



Cl4 \ Xa + ^OXa 



(25) 



where 



Xa=Va + -eVa. 



(26) 



Since we are interested in vector perturbations, we have 
dropped in the above the scalar perturbation Da9- Mul- 
tiplying gives rise to the background equation which 
we have already derived. Multiplying h^, we obtain 



Cl4 ( Xa + -^Oxa 



- (cisD^^afc + ciD^Va + c^DbBaV'') + -0K,(27) 

where we used Eq. (|16p . At linear order scalar- type 
quantity p and p have the same expression as in the 
background. The heat-flux vector qa and the anisotropic 
stress TTab of the aether are given respectively by 



qa = Cl4 ( Xa + -^OXa 

= -cisD'' [aab + D(„Fb)] , (28) 

nab = Ci3 {(Tab + 0<Tab + [D(„VJ,)] ' + OD^,,Vb) } , (29) 



XVa 



where we used Eq. (jl6p to remove A. One can check that 
the momentum conservation equation (jlSp is automati- 
cally satisfied. 

Combining Eq. (j28p with the momentum constraint 
equation. 



D aab = KQa + Qa, 



we obtain 
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1 + Ci3 



[^qa - Ci3D''D(,14)] 



(30) 



(31) 



Note in passing that D^DaVfc = at linear order. 

To proceed further, it is convenient to introduce the 
transverse eigenfunctions. The definitions and the basic 
properties of the eigenfunctions are presented in (28j . In 
terms of the eigenfunctions the vector perturbations can 
be expanded as 



(32) 



where k is the eigenvalue. In terms of the harmonic co- 
efficients, our perturbation equations are written as 



q' + AUq + -H = 0, 



fc^cr 



(2tiS^q - Ci3fcV) , 



1 + C13 

ci4 [V" + 2'HV' + {n^ + W) V] 

+a {U^ -U')V + cik^V 



(33) 
(34) 



-^fc^a, (35) 
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where q := J2i 1i and 11 := ^i- 

Equation psp shows that the fluctuation of the aether 
obeys the wave equation which is similar to the evolu- 
tion equation for cosmological tensor perturbations. The 
crucial difference is the effective mass term which is de- 
pendent on the expansion rate H and the model param- 
eters. The fluctuation of the aether is related to cr via 
the momentum constraint, which in turn translates to 
the magnetic Weyl tensor Hab via Hab = curl dab, and 
thus produces perturbations of geometry. If we neglect 
the matter contents, the above equations reduce to those 
derived in [l^ . 



The velocity Vi of each fluid component is given by 
= Qi/iP^^^ +p^^^). The baryons are coupled to photons 
via Thomson scattering. The baryon velocity Vb obeys 

v't^ + Tivb ^ - — SueCTT (^Vb~ h] , (36) 
Pb \3 ) 

with /i = 4w-y/3. 

We now replicate the photon multipole equations and 
the polarization multipole equations for vectors presented 
in m, m, HI] . The photon multipole equations for vec- 
tors are 



where C := 8/2/4 — and I2 = Il-y/p^, while the polarization multipole equations for vectors are 



B 



e+i 



-kEj 



2£ + r^'-^ 

£ 



(€-Hl)3(2^-M) 2£ + l 



r 



£{£ + !) 



-kB^ 



-SrieCTT 
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£2 



£{£ + !) 



kEf = —Sn^aTB, , 



(37) 

(38) 
(39) 



where Eg and B( are moments of the E and B polar- 
ization.^ The integral solutions to these equations arc 
given in [2^. We will only use the solution for Bi in the 
following discussion: 



S,(7?o) 



1 



dr;fe-^«',[fc(77o-?7)]C, (40) 



where ^'^(a;) := £ji{x)/x corresponding to the projec- 
tion function /3^"'^' in the total angular momentum ap- 
proach [2J], and r represents the optical depth: r := 
J''" ArjSn^aT- The neutrino multipole equations are of 
the form p7[) without the Thomson scattering terms. 



IV. INITIAL CONDITIONS 

We solve the relevant set of equations at early times 
in order to clarify the initial conditions. This is done 
by a series expansion in terms of the conformal time 77, 
following [2^ but now taking into account the presence 
of the aether. 

First, by invoking the tight coupling approximation for 
the baryons and photons it is easy to obtain 



Vb 



Wo 



1 + 



(41) 



^ Here we have corrected the typo found in [28 



where vq is the initial value and R = Spb/ipy. 

The Friedmann equation gives the scale factor in terms 
of the conformal time as 



1 



C ( ,-^22, 



(42) 



where flu := fl^ + £7^ and uj := flm'H-Q / y/Tlj^ . The defi- 
nitions of are the standard ones, namely, densities in 
units of the critical density. 

Neglecting the 0{k^) terms at early times in the 
radiation-dominated era, the perturbed equation of mo- 
tion for the aether ((35l) can be solved to give 



V = Akv" 



1 



1 



v\ iurj 



+^i:V'-'^ + -- - , (43) 



where Ak and -41'^' may depend on k, and 



-1 + ^1 - 8a/ci4 



(44) 



As we are interested in a non-decaying regular mode, 



we assume that a/cn < and set A 



0. In order 



for scalar isocurvature modes not to grow, the condition 
a/ci4 > —1 must be imposed [l^. We advocate this 
constraint and consider the range 



< < 1. 
The neutrino multipole equations read 



0, /. 



2k 



15 



-k(7. 



(45) 



(46) 
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4w^/3 and 



^1^/ Pv Wc use the mul- 



where / 

tipole equations (|46p and the momentum constraint p4p 
to get the foUowing early time solution: 



15 Lorj 

Yarz + 15 

l^* Ci3 



I/* + 4i?* 1 + ci3 



6fe 4i?* 



4 i?; 



1 



3R* 



3 RZ\ ' 15 + 4i?; 



Cl3 



15(1 + I/) i^* + 4i?* 1 + ci3 
In the above we defined 



(47) 
(48) 
(49) 

Akk7]'+''. (50) 



R* := 



z.*:=|(l + i.)(2 + :.), (51) 



where i?i's are defined in the same way as in [28|: Ri, := 
^u/^R, R-y ~ ^-y/^Rj and Rb = flb/flm- The mode 
associated with Bk is identified as the regular vector mode 
in the presence of the neutrino anisotropic stress [28| . 
Since Bk may be fixed independently of the effect of the 
aether Ak, we discard this mode for clarity and focus on 
the initial condition with Ak and Bk = 0. 

Once the inflation model and the subsequent reheating 
history are specified, one can determine the primordial 
spectrum of the vector perturbation and hence Ak ■ Dur- 
ing infiation with e := 1 — H' /H^ = const, one finds, on 
superhorizon scales, that [l5| 



V ' 



1 i-v) 



PI 



(n„-3)/2 



where 



n„ ;= 3 - Jl - 



4e 



(1 - 



(52) 



(53) 



For -1 < a/ci4 < 0, we have 2 - 2e/(l - e) < n„ < 2. 
At the end of inflation, r] = rjg, we have the estimate 



H 

w 



PI 



„/2 



(54) 



where H is the inflationary Hubble scale and corre- 
sponds to the horizon scale at 7] ~ rjc. The factor fc""/^ 
reflects the fact that V decays during inflation if n > 0. 
The amplitude may further change from ([5^ during the 
reheating stage, and hence the primordial amplitude de- 
pends also on the detailed history of reheating. In our 
actual calculation, wc simply assume that 



(55) 



y 0.0 



+ 

^ 0.001 




Tensor TT (r=0.1 
Tensor BB(r=0.1 
Vector TT(nv=2) 
Vector BB (n„=2) 



Multipole 1 

FIG. 1: CMB B-mode polarization and temperature 
anisotropy power spectra in the EA theory. For compari- 
son, those from the tensor perturbation in standard GR are 
also plotted in the case of the tensor-to-scalar ratio r = 0.1. 
In this figure, ci = — 0.2, cia — — 0.3, Ci4 = —a — —0.2, and 
dimensionless primordial power spectra are Vv oc A;"^ and 
Pt oc k°. 



where Ao is a constant. 

As was already derived in Eq. (|43| . V grows on su- 
perhorizon scales, V ~ rj'^ , in the radiation-dominated 
stage. In the matter-dominated stage it turns out that 
the superhorizon behavior is given by I^ '-^ ry"™ with 



-3+ v/l~~24^Vci4 



(56) 



Since — 1 < t'm < 1, 1^ may grow or decay in the matter- 
dominated stage, depending on a/ci4. Therefore, al- 
though the amplitude ([54|) (or Ao) may be tiny as a con- 
sequence of the decay of V during inflation, in the subse- 
quent radiation- and matter-dominated stages V can be 
amplified on superhorizon scales, leading to an obscrva- 
tionally relevant aether perturbation. The special case 
a = — ci4 is an example of such situations, for which 

V = I^ni = 1. 



V. NUMERICAL RESULTS 

Using all the ingredients derived above and the C AMB 
code [33 modified so as to incorporate the presence of the 
aether, we have completed the numerical calculation for 
the B-modc polarization power spectra in the EA theory. 
An example of our numerical results is presented in the 
Fig. [TJ For comparison, we show contributions from in- 
flationary gravitational waves in GR, assuming that the 
tensor-to-scalar ratio is given by r = 0.1. In fact, the 
aether modifies the behavior of the tensor modes as well. 
However, the effects of the aether on the tensor modes 
arc just to shift the location of the peak and to change 
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Tensor TT{r=0.1) 
Tensor BB(r=0.1) 
Vector TT{n„=2) 
Vector BB (nv=2) 



"V"- 




Multipole 1 

FIG. 2: CMB B-mode polarization and temperature 
anisotropy power spectra in the EA theory. For comparison, 
those from the tensor perturbation in standard GR are also 
plotted in the case of the tensor-to-scalar ratio r = 0.1. In 
this figure, ci = -0.019, C13 = -0.03, C14 = -q = -0.0128, 
and dimensionless primordial power spectra are Vv oc fc"" 
and Tt oc kP . 



the absolute amplitude of the primordial spectrum, and 
they are very small for small values of c^. The ampli- 
tude Aq is adjusted so that the low-^ TT spectrum from 
the vector perturbation has the same magnitude as this 
primordial tensor contribution. We see in this case that 
the BB spectrum in the EA theory is larger than that 
from primordial tensor modes at ^ > 100, and hence the 
B-mode is potentially detectable in future CMB obser- 
vations aiming to detect r = O(O.l) - 0(0.01) [H]. In 
plotting Fig. [1] we chose a = — C14. In this case, it has 
been discussed in [isj that the TT power spectrum has 
roughly the same shape as the one from the inflationary 
gravitational waves. As one can sec. Fig. [T] shows the 
same scalings for the two TT spectra. 

Using a more realistic parameter set evading all the 
existing constraints (see Appendix A), we plot the B- 
mode spectrum in Fig. [21 From this we conclude that, for 
a viable range of the model parameters, the B-mode from 
the vector perturbation is potentially detectable in future 
CMB probes even if its amplitude at the end of inflation 
is very small. Note that the amplification of V after 
inflation is determined basically by the ratio between a 
and C14. This means that, even if the model parameters 
Ci are too small to discriminate the EA theory from GR 
with the other observational and experimental tests, the 
CMB B-mode polarization could be a powerful probe for 
the aether field. 

We show the evolution of each variable in a normal plot 
(Fig. [3]) and in a log plot (Fig. |4]). From these figures, 
we find that compared with the aether perturbation V 
and the shear a, matter components are negligibly small 
especially at early times. It can be seen from Fig. 0] 
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log (Scale Factor) 

FIG. 3: Evolution of each variable in a normal plot. The 
horizontal axis is the scale factor in log plot and ao — 1. In 
this figure, ci = —0.2, cm = —0.3, Ci4 — —a = —0.2. 
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FIG. 4: Evolution of each variable in a log plot. The hori- 
zontal axis is the scale factor in log plot and ao = 1. In this 
figure, ci = —0.2, C13 = —0.3, C14 — —a — —0.2. 



that the growth rate of V on superhorizon scales is given 
by y oc S" for a = — C14. This confirms the early time 
solution derived in the previous section. 



VI. ANALYTIC ESTIMATES 

In this section, let us try to understand the shape of the 
B-mode angular power spectrum Cf^ in the EA theory 
in an analytic way. The following discussion is similar to 
the one introduced in [s^ for the B-mode spectrum from 
the inflationary gravitational wave. 

The starting point is the integral solution for the Bg 
which was alreadly introduced in Eq. (|40|) : 



1 



e + 1 



drjie '^'^i[k{ria - ?7)]C- 



Using the approximation for the visibility function. 
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re""^ 5[rj — r]fi), wc see that it is important to know <^ 
at the last scattering surface for determining the B-mode 
polarization. 

Now, wc expand the multipole equations (|37l) , ([38]) and 
([55)1 in terms of fc/f. Neglecting Iq, Ig {i > 3), and Eg 
{£ 2), we have 



2^2 = --Ul-g-'- 



2k 



—C H ka, 

15V 15 



C 

15^ 



At leading order, Eq. (|59p reduces to 



15 



C ^ £^2 



16 



:/2 



(57) 
(58) 
(59) 

(60) 



subhorizon scales {cykr] ^ 1) V decays a.s S ^ with os- 
cillations cos(fcc^77). 

Bearing these facts in mind, we can derive the 
wavenumber dependence of V at recombination r]rcc- Su- 
perhorizon modes evolve in the same way and keep their 
dependence on k, Ak, until horizon crossing. Therefore, 
the superhorizon modes at rj^ec < l/c^r^roc) keep the 
primordial fc-dependence, Ak- After horizon crossing, V 
decays as so that V{riycc)/V{ii^,) oc S'(77*), where 

77* := l/cyk. The modes with l/ci,rircc < k < l/c„77eq 
reenter the horizon in the matter-dominated stage, where 
Tjcq rcferes to the radiation-matter equality time. For 
these modes, we have S{r]^,) oc cx k~^ and 1^(?/*) c>c 
?7^" cx fc""", and hence V{rii-cc) oc k~'^'''''^Ak- Similarly, 
the modes with l/cuTyoq < k reenter the horizon in the 
radiation-dominated stage and so ^(Tyroc) k~^~'^ Ak- 
in summary, the wavenumber dependence of V{ri,-cc) 



Substituting this into Eq. ([55)) and neglecting a higher 
order term, we find 



and 



32 A: / 4 



^ 4 2 5 7- ' 



—a 
3 



(61) 



(62) 



The initial conditions we have adopted in Sec. IV imply 
that at the early times Ji = 41)^/3 = q-y/ Pj and Vi, can be 
neglected compared with a. Then, we have the relation: 



4 k 
15 T 



(63) 



Under the same approximation, a and the fluctuation in 
the aether field V are related as 



Ak (k < l/CyTj^cc), 

V{'qrcc) oc { k^^-^'^Ak (1/Ct,77rcc < fc < l/c„77oq), 
k-^-^Ak (1/C„77cq < k)- 

(67) 

In the above argument we ignored the subhorizon oscil- 
lation of V- 

In all the above calculations, we have neglected the 
matter velocities Vi compared with a. Our numerical re- 
sults in Fig. [3] justify the approximation. The amplitudes 
of V and a are large enough to neglect Vi long before the 
recombination epoch. In Fig. SI we can see that the evo- 
lution of (7 tracks that of V and the ratio is approximately 
equal to —013/(1 + C13) until the neutrino velocity v^, be- 
comes comparable. (In Fig. |4l C13 = —0.3, so that the 
ratio is 3/7.) 

Wc are now in position to derive the simple scaling 
relation for Cf^ in an analytic way. The CMB B-mode 
power spectrum is roughly expressed as 



Cl3 



1 + Ci3 



-V- 



(64) 



The right-hand side of the evolution equation for the 
aether field ([35]) is expressed by the aether field itself 
and this equation becomes the closed form. 



BB 



dlnkrv{k)BiBi, 



(68) 



where Vv{k) is the dimensionless primordial power spec- 
trum for V- From the above discussion, we see 



V" 



2nv' + clk^v 



1 



Cl4 



1 



a 

Cl4 



where the sound speed is given by 

^2 



.2 _ £L 

Ci4 



1 - 



'13 



n' 



2ci(l-Hci3) 



y = 0, (65) 



(66) 



Equation ([65]) tells us that on superhorizon scales 
[cykr] ^ 1) F is proportional to rf (or, equivalently, S^) 
in the radiation-dominated stage and to rf^ (or, equiv- 
alently, 5'^"/^) in the matter-dominated stage, while on 



'°AXf , ije[k{vo - v)] k V ci3 

kivo -V) T Akl + ci3 

ije[kiVO - '7rec)] k V{rirec) Ci3 



k{T]o - ?7roc) T Ak 1 + Ci3 ■ 



(69) 



Using the fact that the projection factor £je{x)/x has a 
peak at £ '-^ x, the power spectrum reduces to 



^kj k= — ^ 

■no -Vrcc 



fc2+"^[M',(fc(77o-7?rcc))]'dlllfc. 

(70) 

For example, taking n„ — 1, the above integral can be 
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evaluated to give the scaling 



0\-1v 



(4cc < I < 4q), (71) 
(£eq < ^ < ^A^„J, 

showing a peak at £poak ~ ^roc, where we have defined 
4oc := (??o - rjiec) / CvVrcc and 4q := (vo - V^oc) / Cyrjcq- 
^Ariioa represents the scale over which the phase-damping 
effect shows up due to the width of the last-scattering 
surface. On scales £ > ^Ar^oci ^^'^ above approximation 
is no longer justified. In deriving the scaling we used the 
following integral formula for the Bessel function, 



dlnfcfc2+™*i(fc) 
,2v^ r(l-m/2) 



T{l + m/2) 



4 r(3/2-m/2)r(; + 2-m/2)' 



(72) 



The scaling behavior can indeed be seen in Fig. [5] for 
the illustrative case a = — C14 {v ~ i^m ~ 1), though the 
scaling in the range £cq < £ < i /\n^ is not clearly seen 
due to the phase-damping effect [39[. Actually, in order 
to evaluate the shape at the smallest scales correctly, we 
must take into account more complicated physics such as 
neutrino anisotropic stresses and the Silk-damping effect. 
We would emphasize that our numerical calculations in- 
corporate all of these effects. For instance, in Fig. [5] we 
can confirm the oscillatory behavior at £ > ipcak arising 
from the subhorizon oscillation of V. (In Fig [5l the B- 
mode spectrum from tensor perturbations in GR is also 
plotted for comparison.) 

We can also gain an understanding of how the shape 
of the angular power spectrum depends on the model 
parameters. From Fig. [SI one can confirm the follow- 
ing three things: (i) Since the angular power spectrum 
on the largest scales £ < £rcc depends only on the pri- 
mordial spectrum, the plotted examples show the same 
scaling at this scale. This means that the observation of 
the polarization at the largest scale can determine the 
primordial spectral index; (ii) The peak position is in- 
versely proportional to the sound velocity of the aether 
vector perturbation Cy. Actually, if we find the peak of 
the observed B-mode spectrum, the information is di- 
rectly converted into the exact value of Cy in case the 
other cosmological parameters are already determined; 
(iii) The difference of the small scale scaling arises due 
to the difference of the growth rate of V on supcrhorizon 
scales, which can be seen in Eq. ([7T|) . 



VII. SUMMARY 

In this paper, we have considered the Einstein-aether 
theory in which gravity is modified through the addi- 
tional vector degree of freedom (the aether), and studied 
possible signatures of the aether field in the CMB polar- 
ization. In standard GR, vector cosmological perturba- 
tions simply decay without sources, and hence they are 



B3 



Vector BB(n„=1) 




TensorBB (r=0.1) 
















..■'/ / -^peak 









Multipole 1 

FIG. 5: (a) Scaling for the illustrative case with n„ = 1 
and a = — C14 = 0.2. The other parameters are given by 
C13 = —0.3, and ci = —0.1; (b) Parameter dependence of the 
spectrum. In the two examples C14 is different while the other 
parameters are fixed as n„ = 1, C13 = —0.3, and c\ = —0.1. 
The primordial amplitudes are arbitrary. 
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FIG. 6: (a) Scaling for the illustrative case with — 1 
and a — — C14 = 0.2. The other parameters are given by 
C13 — —0.3, and ci = —0.1; (b) Parameter dependence of the 
spectrum. In the two examples C14 is different while the other 
parameters are fixed as n„ = 1, C13 = —0.3, and ci = —0.1. 
The primordial amplitudes are arbitrary. 



less relevant to observations. In the Einstein- aether the- 
ory, however, the vector modes are dynamical, so that 
they could leave imprints on the CMB signals in a way 
similar to the infiationary gravitational waves. 

Using the CAMB code modified to incorporate the 
aether vector perturbation, we have computed the CMB 
B-mode polarization power spectrum. We have found 
that the amplitude of the B-mode polarization from the 
aether vector can be larger than that from inflationary 
gravitational waves with r = 0(0.1) on small angular 
scales, which would be measurable with near future CMB 
probes. 

Moreover, we have found that, for a set of parameters 
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Ci evading all the existing observational and experimen- 
tal constraints, the B-mode polarization spectrum indeed 
shows distinguishable features from that from inflation- 
ary gravitational waves in GR. Thus, the B-mode polar- 
ization spectrum would potentially be a crucial test for 
the Einstein-aether theory. 

We performed analytical calculations, by which we 
clarified the physical origins of the shape of the B-mode 
power spectrum. The scaling relations derived in an ana- 
lytical way turned out to be consistent with the numerical 
results. 

In this paper, we have simply assumed the primor- 
dial amplitude of the vector perturbation which is gen- 
erated anyway in the presence of the aether field. The 
primordial spectrum actually depends upon the underly- 
ing inflation model and the reheating process. It would 
be interesting to explore whether the B-mode signal is 
measurable or not assuming concrete infiation models. 
We will report elsewhere the evolution of vector pertur- 
bation during the whole history of the Universe starting 
from inflation in the EA theory. 



The easiest way to pass the stringent observational con- 
straints, Qfi < 10~^ and 02 4 X 10~^ [i^, is to set ai,2 
exactly to zero by imposing the conditions 



C2 = 



C4 



-2cl - C1C3 
3ci 



Cl 



(A3) 
(A4) 



which is possible since EA theory has four free parame- 
ters Ci{i = 1 — 4). Under these conditions, all the PPN 
parameters in the EA theory coincide with those of GR. 



2. Stability of each perturbation mode 

Linear perturbations around the flat or FRW metric 
have been studied in [isl, H^l • They concluded that quan- 
tum and classical stabilities of tensors (spin-2 modes), 
vectors (spin-1), and scalars (spin-0) constrain the aether 
parameters to be 
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Cl3 > -1, 
-2 < Ci4 < 0, C123 < 0, 

2ci < c?3(l -hcia). 



(A5) 
(A6) 
(A7) 



3. Radiation damping and strong self-field effects 



Appendix A: Observational constraint 

In this appendix, we summarize the existing observa- 
tional constraints on the aether parameters ci , • • • ,04, 
following ^\^. 



1. Post-Newtonian limits 

Parametrized post-Newtonian (PPN) parameters in 
the EA theory have already been analyzed in [4l|, |4^ . 
Two PPN parameters, the Eddington-Robertson-Schiff 
parameters /? and 7, are identical to those in pure 
GR [4l|. The Whitehead parameter, ^, which charac- 
terizes a peculiar sort of three-body interaction vanishes 
in the EA theory [l^], and the five energy- momentum 
conservation parameters and Ci,2,3,4 vanish because 
the theory is derived from a Lagrangian. 

The aether defines a preferred frame, and its effect is 
encoded in the remained PPN parameters ai and a2- 
The exact values of ai and a2 arc found in [42I : 



ai 



-8(ci 



C1C4) 



C(2 — -Z- 



2ci 
2 



(2ci3 - ci4){a + C14) 

Ci23(2 - C14) 



(AI) 
(A2) 



The radiation damping rate in the weak field limit was 
first calculated in [43, and then strong field effects are 
included in [4^ . They found that in general the damping 
rate in the EA theory is different from that of GR. How- 
ever, in the special case with ai = 02 = 0, it is identical 
provided that the quadrupole coefficient of the radiation 
damping rate. 



A 



Cl4\ 

2 / 



Ci4 



6(2-hci4) 



1 

Ct 



2ci4c2 



13 



1 



(2C1 -|-Cl3C_)2 



2a2 



ai 



2(2ci3 - C14) 



1 



(A8) 



is equal to one. Here we defined the new parameters 



1 



1 + Ci3 
2 (2 + Ci4)ci23 



(l + ci3)(2-a)ci4' 



(A9) 
(AID) 



which correspond to the velocity of tensor modes and 
scalar mode, respectively, and c_ := ci — C3. 

Foster [i^ derived several constraints by considering 
strong field effects. All of them reduce to the condition 
\ci\ < 0(0.1) with = a2 = under the current obser- 
vational uncertainties. 
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4. Cherenkov Radiation 

If the sound speed of each perturbation mode were 
smaher than the speed of Hght, then the ultra-high- 
energy particles would loose their energies into the mode 
in a similar way to Cherenkov radiation. Using this fact, 
Elliott et al. [43 derived limits on the aether parameters. 
If all the aether modes propagate superluminally, wc do 
not need to take into account their constraints. 

The conditions for superluminal propagation of ten- 
sors, vectors, and scalars are 



C13 < 0, 

(2 -I- Ci4)ci23 < (2 - a)(l + Ci3)ci4, 



2C4 > 



-13 



/(1 + C13), 



(All) 
(A12) 
(A13) 



respectively. The connection between superluminality 
and violation of causality has been under debates |48l - 
[52j . Since the above conditions together with the PPN 
vanishing conditions (|A3|) and (|A4p are equivalent to the 
stability condition of scalar and vector modes, we here 
allow for the superluminal propagation in the EA theory. 



5. Scalar mode constraint 

Armendariz-Picon et al. [lH| have formulated the evo- 
lution equations for perturbations around the FRW 
background metric and calculated CMB temperature 
anisotropy spectrum. They found mainly two constraints 
on the aether parameters: 

• In order for the scalar isocurvature mode not to 
grow on superhorizon scales. 



• In order not to have too large an anisotropic stress, 
|ci3| < 1 (A15) 

must be satisfied. 

Appendix B: Allowed parameter values 

If one imposes that the PPN parameters coincide ex- 
actly with those in GR, one has two constraints derived 
from ai ~ a2 = 0. We are now left with two free param- 
eters, for which we use C13 = ci + C3 and c_ = ci — C3. 
Interestingly, the special combination a = — C14 is auto- 
matically satisfied, and we have 



a = — Ci4 = — 2- 



C13C- 



(Bl) 



Ci3 -h C_ 

Imposing that the perturbations are stable and super- 
luminal, we can restrict the two parameters within the 
range as 



Cl3 



3(l + ci3) 



1< Ci3 <0, 
< c- < 0. 



(B2) 
(B3) 



a < — ci4 



(A14) 



For the parameters satisfying the above constraints, we 
have safely a cosmological solution (see Eq. (HU), and 
we we do not have a growing isocurvature mode (see Eq. 

mi) ). 

Imposing further that ^ = 1, one can evade the con- 
straint from radiation damping rate. Thus, we are finally 
left with a single parameter, say, C13, within the range 
dm 0. Taking \ci\ < 0(0.1), we choose to use 

ci = -0.019, C2 = 0.014, C3 = -0.011, C4 = 0.0063, 
Cy = 1.241, a = -ci4 = 0.0128 (B4) 



must be satisfied. 



in Sec. V of the main text. 
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